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, Abstract 

■ In a series of papers we have developed the method of iterated sector decomposition 

\ for the calculation of infrared divergent multi-loop integrals. Here we apply it to phase 

• space integrals to calculate a contribution to the double real emission part of the e+e^ — > 

I \ 2 jets cross section at NNLO. The explicit cancellation of infrared poles upon summation 
over all possible cuts of a given topology is worked out in detail for a specific example. 
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1 Introduction 



The precision measurements at LEP, SLAG, HERA and the Tevatron in the past years made 
it obvious that QCD corrections at next-to-leading order (NLO) accuracy are mandatory for 
a successful comparison of data and theory. Fortunately, a multitude of NLO corrections 
are available these days and the techniques for these calculations have reached a rather 
mature state, at least in what concerns 1 ^ 3 or 2 ^ 2 partonic processes. 

However, the advent of the LHC and above all an e^e~ linear collider will provide 
instances where the NLO QCD corrections are - in some particular cases - not sufficiently 
accurate to match the experimental precision. This is true for example for the extraction 
of the strong coupling constant ag from jet observables at a future e^e~ collider. 

This fact served as a motivation for substantial progress towards the calculation of 
NNLO corrections to important processes, in particular in what concerns the two-loop 
virtual corrections CP-IH] and the one-loop corrections combined with real radiation where 
one additional parton is emitted [3 IHl IHl- The construction of a fully differential NNLO 
partonic Monte Carlo program however also requires the calculation of the real emission part 
where up to two partons can be unresolved, and - last but not least - the combination of all 
the ingredients to a stable and sufficiently fast Monte Carlo program. In what concerns the 
double real emission, subtraction schemes have been proposed in the literature |in |llll[T2] . 
but a complete calculation including the final Monte Carlo program has been performed 
so far only for the particular case of the photon + jet- rate in e^e~ annihilation |13| I14j. 
However, a lot of activity concerning this subject is going on at the moment. The efforts are 
concentrated particularly on the process e'^e~ Sjets, as this reaction is both appealing 
from a phenomenological point of view (e.g. measurement of ag) as well as from a theoretical 
one (no problems due to initial state singularities). Nevertheless, it is worthwhile to consider 
first the reaction e~^e~ 2 jets as this is the simplest example where the treatment of the 
double unresolved real radiation can be studied, and therefore a good testing ground for a 
new method. 

Progress in what concerns the integration of subtraction terms has been made in |15j . 
where it has been shown that the integrals of any 1 — > 4 matrix element in massless QCD 
over the total phase space can be expressed by four master integrals. These integrals have 
been evaluated analytically as well as numerically. 

The complexity of the phase space integrals with two unresolved partons stems from the 
fact that the corresponding IR singularity structure is overlapping. We have demonstrated 
in ii-ow one can disentangle in an automated way overlapping singularities in parameter 
representations of dimensionally regulated multi-loop integrals. The method of iterated 
sector decomposition, combined with numerical integration of the pole coefficients, proved 
successful to deal with very complicated Feynman diagrams |161 117j . 

In this paper we apply the same ideas to parameter representations of phase space 
integrals, as a continuation of the work done in jl5l I18j . In the same context the method 
of sector decomposition also has been applied meanwhile in |19j . 

In contrast to jl5j , where this method of numerical integration has been applied to 
master integrals only, we show here that we also can deal with numerators coming from 
gauge couplings and thus finally with complete matrix elements. Working through one 
example in great detail, we demonstrate the viability of our approach, showing explicitly 
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Figure 1: Sample three loop topology related to e"^e 2jets at NNLO. 

the cancellation of infrared poles for a given topology when summing over all cuts. 

The paper is organised as follows. In section 2 we discuss in detail the cancellation 
of infrared divergences in a concrete example. We work out the cut structure and the 
counterterms and represent the considered topology by a number of cut diagrams which 
give rise to phase space integrals. In section 3 we present the analytical results for the 2- 
and 3-particle cut diagrams, and in section 4 we discuss the method to perform the phase 
space integration of the 4-particle cut numerically. In section 5 we present the numerical 
result for a more complicated topology which has up to poles. Section 6 closes the 
paper with a discussion. 

2 Cancellation of infrared divergences 

Two-jet production in e^e~ collisions at NNLO in corresponds diagrammatically to sums 
over cuts of three-loop vacuum polarisation graphs. After renormalization of ultraviolet 
sub divergences the three-loop diagram has at most an overall ultraviolet divergence with a 
finite imaginary part. Because this imaginary part is related to the sum over all cuts of the 
diagram, one concludes that the sum over all cuts has to be finite. This is just the KLN [20] 
cancellation mechanism for infrared final state singularities. In our case, one has 2-, 3- and 
4-particle cuts which correspond to phase space integrations over the respective partons. 
We require 2 jets in the final state, such that up to two partons can become unresolved. 

Consider for example the diagram shown in Fig.^ In the following subsections, we will 
work out the UV renormalization and cut structure of this diagram. We use Feynman gauge 
throughout the paper. 

2.1 Cut diagrams 

Infrared cancellations take place when summing over cuts of a given renormalized topology. 
In our case we have the 2-, 3- and 4-particle cuts 02,0^ and C4, which graphically are 
given by 
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The factors zi and Z2 denote the contributions from counterterms. The sum of terms 
in each Une in the above figures is formally UV finite. In dimensional regularization the 
radiative corrections to the two-point functions vanish if the particle is put on-shell. This 
is the mechanism which formally leads to the conversion of UV into IR poles. Using the 
diagrammatic rule 




the only remaining terms in the sum C = C2 + C3 + C4 are 
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(1) 

The 1/e terms contained in zi and Z2 now represent IR poles. The graphs shown in ((T|) 
denote the contribution from the given topology to the full process. As the imaginary part 
of the corresponding 3-loop topology is finite, unitarity implies that C is IR finite. 

The part which is hard to calculate analytically is the 4-particle cut, that is why we 
advocate a numerical evaluation of C4. The remaining combination of counterterms and 
2- and 3-particle phase space integrations can straightforwardly be done analytically. This 
will be worked out in the next section. 



2.2 UV renormalization 

To compute the renormalization constants zi and Z2 we have to determine the pole parts of 
all graphs and subgraphs. Following a graphical BPHZ notation one has to carry out the 
following subtractions to renormalize the one-loop selfenergy and finally the whole diagram: 
The one-loop subtraction: 




= finite 



and the two-loop subtractions: 
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finite 



The shaded boxes denote the MS prescription to keep only the pole part of the expres- 
sion, up to absorption of a standard factor into the coupling: 
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Using the integral 

I{a,(3) 



iT,D/2 [_fc2]a[_(^_p)2]/3 

^(2-Z))(V) 



2^D/2-»-g r(a + P- D/2) T{D/2 - a + l)T{D/2 - (3) 
T{a)T{(3) T{D -a- (3 + 1) 



(3) 



one finds the following analytic results for the graphical expressions [D = A — 2e). At one 
loop: 



a. 



r(l - e)2 (1 - e) 



\i7rj \ /i^ / r(2 - 2e) e 



and for the two-loop part: 
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The renormalization constants can be read off directly: 



zi = Cp 
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(4) 



(5) 



(6) 



(7) 



The non-local logarithmic terms cancel, as guaranteed by the BPHZ theorem. 

3 Matrix elements and phase space integrals 

What remains to be done is the evaluation of the phase space integrals for the cuts shown 
in Eq. The corresponding matrix elements are given by the following formulae, where 
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Pi (P2) are the momenta of the quark (anti-quark), ^3 and p4 denote the gluons. 

l^'-l' = «(^')'(i-')'"°.c^( r(i;'K4,J S 

/ 2e \ 3 



^r(l + e)(47r) 

(Sl3 + S34)(S12 + ■§24) - S23 



(10) 



|7Wi^2p and |A^i^3p can be integrated directly. The corresponding phase space integrals 
are given in detail in the Appendix. We obtain 

T,_. ^ /.*._.|A.._,|= ^ 2.0= f^i^l^ (U) 

Using the integral 

00 2 

M^) = f dyidy2dy3 5{l-yi-y2-y3){yiy2y3)''— = -— — — (12) 
J y2 e r(3 - 3e) 

where yi = Si2/Q'^,y2 = sis/Q^iys = ■523/Q^, one finds for the 1^3 case: 

|2 

-2e 

87r3 ^' I r(l + e)2 



Ti^3 = J d^i^3\Mi^3\ 



r^r. r r^'Y ^(3 - 26) X^(2 - 26) 



_ /Qn- 2(i-6)2r(i-6)2 
- -"^^^-^l^^j r(i + 6)r(3-36) 

The 1 — > 4 case cannot be solved easily analytically. We write it in the following form: 

J d$i^4|A^1^4p 



■1^4 



where the nontrivial integral J^ie) remains, which will be dealt with numerically in the next 
section. 

So far, one finds for the cut diagrams in Eq. (^: 

2~ n 1 



C2 = ztT,^2[-^--^e 
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(15) 



and consequently the poles contained in J4(e) have to cancel with 



C2 + C3 I 

pole part — T\^2 
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4 Numerical evaluation of 1 ^ 4 phase space integrals 

In massless QCD, the integrals of any 1^4 matrix element over the total phase space can be 
expressed by four master integrals whose analytical evaluation is complicated, but has been 
achieved in ^5]. This allows to follow the conventional procedure to establish a subtraction 
scheme and to integrate the subtraction terms analytically. However, as the finite part of the 
phase space will finally be integrated numerically anyway, a flexible, completely numerical 
method would be welcome. Of course the problem consists in the isolation and subtraction 
of the infrared poles, stemming from the integration over unresolved particles, before a 
numerical evaluation is possible. This is where sector decomposition is very convenient. 
How it proceeds is shown for the sample integral J4 defined in eq. (|17() . Note that this 
integral contains an integrable singularity of square-root type. For the sake of numerical 
stability it is preferable to perform a mapping such that the integrand is bounded near the 
phase space limits. In detail we proceed as follows: 

Our starting point is the integral J4(e) which represents (up to an overall factor, see 
eqs. (fTn)) . (fT^ ) the integral over the 4~particle cut of the topology in Fig. ^ 



where we have rescaled the Mandelstam invariants by 

and A = \{yiye, 2/2^51 2/32/4) is the Kallen function A(x, y, z) = + y'^ + — 2xy — 2yz — 2xz. 
The derivation of the phase space integral can be found in the appendix, see eq. (|57|l. 

Primary sector decomposition 

To eliminate the delta distribution in (|17j) we split the integration region into 6 "primary 
sectors" by the following decomposition of unity: 



Me) = ^ fUdvi ei-x)i-x)-'/'-^6ii - y: yj) 




(2/1 + 2/5) (2/2 + 2/6) -2/3 2/4 
2/2 (2/2 + 2/4 + 2/6)^ 



(17) 
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(18) 



j=ij^k=i 



7 



In each primary sector j we apply the mapping 




tkVj if j 
Vj if k = j 



(19) 



Note that to each index there exists a conjugate index defined by the pairing of Mandelstam 
variables in the arguments of the Kallen function, X{yiyQ, 2/22/5; 2/32/4)- The constraint 0(— A) 
is solved for the variable with index conjugate to j in each primary sector j. To be explicit 
we give the result for primary sector 1 in the following, where we have 

A(t6; ^2*55 ^3*4) = 4^ tf = t2t5 + t3t4 ± 2^/t2t^t4t5 = (^7*2^5 ± Vhti)"^ 

Remapping the square-root singularity 

Before iterated sector decomposition can be applied to disentangle the IR singularities, it is 
necessary to perform some variable transformations such that finally all possible singularities 
are at zero. Further it is useful to make a quadratic transformation tj = xj for j = 2, 3, 4, 5 
to avoid square roots. 

One possibility of remapping is to split each primary sector into two subsectors A and 
B, with Iq G [ig'jie] = ^ te £ [^e'^e"] = -^i where tg = t2t5 + ^3^4. The following trans- 
formations then lead to a numerically stable behaviour near the phase-space boundaries: 



Note that the 0-function constraint is trivially fulfilled whenever one of the variables Xj, 
j £ {2, 3, 4, 5} goes to zero. However, at this point we have not mapped all possible 
singularities to zero yet. For example, can vanish if {x2, xs, X4, X5} — > 1,X6 — > 0. In this 
case another transformation Xj — > 1 — xj is made^ . Note that all these transformations are 
done automatically in our program. Having finally mapped all possible singularities to zero, 
iterated sector decomposition can be applied straightforwardly. One obtains the coefficients 
of the 1/e poles as finite parameter integrals which can be integrated numerically, in the 
same way as has been explained in jl5l I16j . 

^It can occur that Xj is singular as well for some j € {2, 3, 4, 5}. In this case the integration region for 
Xj is split at 1/2, and only after this splitting the variables are remapped such that all possible singularities 
are located at zero. 




(20) 



(21) 



The function F depends on the topology, in our case 
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Figure 2: Cut topology with poles. 



Numerical solution 

For the integral J4 we find numerically 

(1 - e)2 J4 = o.50o4 + 2.749- + 7.869 + 0(e) (23) 

e 

Inserting this result into (|15j) and comparing to eq. (|16j) we see that the poles are cancelled 
within the numerical precision which we chose to be 0.1%. 



5 Another Topology 

To show that our method also works in the case of a more complicated pole structure and a 
lengthier numerator, we consider now the cut graph shown in Fig. [21 which leads to poles up 
to and the expanded numerator contains about 70 terms. The corresponding matrix 
element is given by {pij... = i>i+ pj + ■ ■ ■) 

,2 n3 2 ^2 ( ^ ^ tr(pi7>i37Ppi347^p27pP247<TP2347M) 

A1i_»4 = — (47r) aa„Lp —— , , , , 

' ' V ; . ^ \^r(l + e)(47r)^y S134S13S234S24 

(24) 

Writing again the phase space integral in the form 



-2t 

IJ? ) 16(l-e)r(l + e)2r(l-2e) 



we obtain numerically 



~ 1.000 3.000 3.776 8.957 „^ , , 

J4 = — ^ + + + — + 57.85 + ©(e) (25) 



We stress that no analytical manipulations whatsoever, in particular no reduction to master 
integrals, are necessary to achieve this result. 



6 Discussion and Outlook 

In this paper we have argued that the method of sector decomposition as proposed in [THIITK] 
and applied in a similar form in FT^ can serve to calculate the double real emission part 



needed for NNLO corrections to jet cross sections in massless QCD. We have shown expHcitly 
how the IR poles cancel by considering one sample topology contributing to e+e" — > 2 jets 
at NNLO. We calculated all cuts contributing to a 2 jet final state: The double real emission 
part, corresponding to a 4-particle cut, is calculated numerically by sector decomposition, 
and all the poles are shown to cancel with contributions from the 2-and 3-particle cuts 
within the chosen numerical precision^. We also calculated numerically the 4-particle cut 
of a topology which has the maximal possible number of IR poles. 
The method has several appealing features: 

• The (overlapping) soft /collinear poles are extracted without the need to establish a 
subtraction scheme and to integrate analytically over complicated subtraction terms. 

• The finite parts are available as regular functions of the kinematic invariants. If a mea- 
surement function is included, as already has been done in one can obtain fully 
differential cross sections, such that for example the implementation of experimental 
cuts should be straightforward. 

• It is justified to expect that the Monte Carlo integration to obtain the final cross sec- 
tion does not lead to major instabilities as the singularities at phase space boundaries 
already have been remapped. 

• The generalisation of the method to the 1 — > 5 phase space is feasible but has to be 
investigated further. 

The complete double real emission part contributing to e~^e~ — > 2 jets can be evaluated 
along the same lines and will be given elsewhere. 
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Appendix 

The D-dimensional 1 ^ phase space for = 2,3,4 

The differential form of the phase space for a 1 — > particle phase space in D dimensions 
is given by 

N N 

d<^i^N = (27r)^-^(^-i) [n d%S{p''j)e{Ej)]6[Q - ^^p,) (26) 

i=i i=i 
Here Q is the incoming momentum and the pj denote outgoing particles with light-like 
momenta and energy component Ej. As Q is time-like for physical kinematics one can 
always achieve Q = {E,0^^~^^) by an adequate Lorentz boost. Let us specialise now to the 
cases A^ = 2, 3, 4. 

^Note that the method of sector decomposition could also serve to calculate the 2-and 3-particle cut 
contributions. In this case the algorithm for multi-loop integrals as described in |16l 1171 can be applied. 
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Case 1 2: 

For N = 2 the momenta can be parametrized by 

Q = (i?,0(^-i)) , pi=i?i (1,0(^-2), 1) , p2=Q-pi (27) 

Integrating out the (5-distributions leads to 

d^i^2 = (2vr)2-^ 2'-^ (Q2)^/2-2 (28) 

where dflD-2 is the differential surface element of the Sd-2 sphere. Its integral is equal to 
the volume of the {D — l)-dimensional unit sphere 

D-l 

dnD~2 = V{D - 1) = (29) 

JD-2 

Case 1-^3: 

For = 3 one can choose a coordinate frame such that 

Q = (S,0(^-l)) 

pi = i?i(i, 0(^-2), 1) 

P2 = S2(l,0^^~^\ sine, COS 0) 

P3 = Q -P2-Pi (30) 
Integrating out the ^-distributions yields 

d$i^3 = ^(27r)3-2D dEidE2dei[EiE2s\iie]^'^dVLD-2 d^o^s (31) 

As in the following a parametrization in terms of the Mandelstam variables Sij = 2pi -pj will 
be useful, we make the transformation Ei,E2,0 — > si2, S23, S13. To work with dimensionless 
variables we define yi = si2/(3^, y2 = ^is/Q^, ys = S23/Q2 which leads to 

d$i^3 = {27rf-'''^{Q')''-'dnn-2dnn-3[yiy2y3]''/'-' 

dyi dy2 dys G(yi) 6(^2) 6(^3) ^(1 - Vi - y2 - Vs) (32) 

Case 1 ^ 4: 

Starting from Eq. (|2())) and eliminating p^ yields 

= (2.)-3z^^^^^e(i^,)e(i^2)e(i^3) 

e(S -E1-E2- Es) 6[iQ -P1-P2- P3?) (33) 
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Choosing a frame where 



Q 






pi 


= (1,0(^-2) 


1) 


P2 


= E,{i,6^^-'^ 


sin 6*1, COS 61) 


P3 


= Es{i,6(^-'^ 


sin ^3 sin 62, cos 


P4 


= Q-P1-P2- 


-P3 ■ 



(34) 



leads to 



1 



[iTif-^^ dEi dE2 dEs dOi d92 de^[EiE2E^ sin 6*1 sin 6*2]^"^ sin (9; 



D-3 , 



3D-4 



df^D-4e(Si) 6(^2) e(^3)e(^ -E1-E2- E3) 

6{E^ - 2E{Ei + E2 + E3) + 2{pi ■P2+P1-P3+P2- P3)) (35) 

As above we map the angle and energy variables to the Mandelstam invariants as integration 
variables. The Jacobian in combination with terms already present in (|35() can be written 
as the determinant of the Gram matrix Gij = 2pi ■ pj. The determinant can be expressed 
by the Kallen function A(x, y, z) = + + — 2xy — 2yz — 2xz as 



det(G) = A(si2S34,Si3S24,Sl4S23) 

= —[4E El E2 E3 sin 9i sin 62 sin 9^ 



(36) 



We see that det(G) has to be negative semi-definite. With the dimensionless variables 

yi = si2/Q^ , y2 = sis/Q^ , ys = S23/Q^ , 2/4 = si4,/Q^ , ys = S24/<3^ , ye = ssa/Q^ 
and A = A(7/i7/6, ^22/5, 2/32/4) one obtains finally 



d^i. 



(2^)4-3D(g2)3D/2-4 2-2D+l^^^_^ 



n dyMyj] 



G(-A)[-A](^-^)/25(l-;^y,) 



(37) 
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